Invariant formulation of the Functional Renormalisation 
Group method for U(n) x U (n) symmetric matrix models 

A. Patkos 
Institute of Physics, Eotvos University 
H-1117, Pazmany Peter setany 1/A, Budapest, Hungary 

<n: October 25, 2012 



(N 

s 

CD 



> 
O 



O 



Abstract 



o 

-t— > 
O 

The Local Potential Approximation (LPA) to the Wetterich-equation is formulated explicitly 
in terms of operators, which are invariant under the U(n) x U(n) symmetry group. Complete 
formulas are presented for the two-flavor (£7(2) x £7(2)) case. The same approach leads to a 
unique natural truncation of the functional driving the renormalisation flow of the potential of 
the three-flavor case (C/(3) x U(3)). The procedure applied to the SU(3) x SU(3) symmetric 
theory, results in an equation, which potentially allows an RG-investigation of the effect of the 
't Hooft term representing the Ua(X) anomaly, disentangled from the other operators. 



1 Introduction and General Analysis 



The Ul(N) x Uf>(n) chiral symmetry is relevant to the QCD of n massless quark flavors. Quark masses 
^t" below Aqcd introduce tolerable explicit symmetry breaking which leaves this symmetry still a valid 
approximate concept underlying the construction of effective low energy models of strong interactions. 
At zero temperature the chiral symmetry is realized in the broken phase. The breakdown originates 
{Nl '. from two sources. It is spontanously broken to Uv{n) [UEIE] an d the axial Ua(X) quantum anomaly 
[4, '5j further breaks it to Uy(l) x SUy(n). Low energy effective meson models are constructed along 
this specific breakdown of chiral symmetry. 

Both the strength of the axial anomaly and the chiral condensate diminish with increasing tem- 
perature and eventually the full Ul{ti) x Ur(ti) symmetry will be restored. It is an interesting 
issue to understand how wide is the domain of attraction of the symmetry breaking pattern ob- 
served in nature. Is there any sign of fine tuning in the observed spectra of light mesons? More 
specifically, one should ask if a slight change of the couplings in the UV (high temperature) regime 
could result in a qualitatively different IR (low temperature) ground state. For instance, for n = 2 
there are two alternative breaking patterns: Ul(2) x Ur(2) — > LV(2) — > Uy(l) x SU(2) and 
Ul(2) x Ur(2) — Y U(l) x U(l) x U(l). Similarly, for n = 3 there exists the possibility for a chain of 
symmetry breaking evolution: Ul(3) x Ur(3) — > Uy{?>) — > Uy{2) x U{1) even if one does not consider 
the effect of the axial anomaly. This alternative scenario was investigated recently with help of a 
large- N approximation [6j [7]. 
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A most adequate tool to answer this kind of questions is provided by the Wilsonian Renormal- 
isation Group analysis [8] of the scale evolution of the effective potential (or free energy at finite 
temperature). Starting with some initial potential at small distances the potential settles in one of 
the minima of the large wavelength effective potential. A global exploration of the dependence of the 
symmetry breaking pattern realised in the absolute minimum on the initial set of couplings decides 
on how robust is the actual scheme. A non-perturbative realisation of this program is offered by the 
various formulations of the Functional Renormalisation Group (FRG) equations [9| \TU \ ITT | [T2]. 

Several studies in the n = 2 case were done concentrating on the symmetry breaking pattern which 
is the closest to what is observed in nature [131 US]. But the full landscape was not explored yet, 
therefore we do not know how generic is this pattern. The effect of the anomaly was not included 
into the investigations, though this becomes particularly interesting in the 77 — 77' phenomenology, 
present only for n — 3. 

There is a large number of FRG-investigations of the two-flavor effective models where the full 
symmetry is restricted to SU{2) x SU{2) ~ 0(4) [161 El H], an d the symmetry breaking pattern 
0(4) —7- 0(3) is investigated also by taking into account confinement effects [12]. In the context of 
cold quantum gases also the 0(4) — > 0(2) x 0(2) breakdown was studied with FRG approach[20j. 

This note gives some technical details necessary to the full exploration of the coupling space for 
the U(n) x U(n) symmetric models, n = 2,3. We present the most convenient form of the RG-flow 
equations which will be formulated explicitly in terms of variables invariant under the symmetry and 
respecting the analiticity features of the potential in these variables. 

In this investigation we use the so-called Local Potential Aproximation of the Wetterich equation 
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with a specific choice of the cut-off function suppressing low momentum modes in the functional 
integration [22J 

R k (q) = (k 2 -q 2 )Q(k 2 -q 2 ). (2) 

It leads to the equation 

h d+l Oj n h d+1 Oj d 

\ ' eigen—modes,a 'a \ I 

where d is the dimension of the space-time, Qd the surface of the d- dimensional unit sphere and 
the sum a runs over the different eigenmodes of the mass-square matrix with eigenvalues M 2 of 
degeneracy g a . The second equality makes it explicit that the scale evolution of the potential is 
driven by the characteristic polynomial n[(A; 2 + M^) 9a ] of the mass-square matrix. 

The potential U k is invariant under the symmetry transformations, therefore it depends only on 
invariant combinations of the field variables. In the course of the derivation of the Wetterich equation 
one assumes that it has a power expansion in terms of these variables. Also the right hand side of 
03]) should share this feature. We refer to this requirement below as the analyticity of the right hand 
side of (J3J. 

The matrix field of the U(n) x U (n) invariant theory is parametrized with help of the generalized 
Gell-Mann matrices as 

$ = 0A', I = s, + i7r, J Z = 0,l,2,...,n 2 -1, Ti\ l \ k = 25 kl . (4) 
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The quadratic and quartic U (n) x U (n) invariants are defined as 



p = Tr$t$ = 2( S f + 7rf), 



1 N 2 

r = Tr $ f $ Tr$ f $ 



n 



(5) 



Usually, one considers only U(p,r), since these are the two operators, perturbatively renormalisable 
around the Gaussian fixed point. The expression of the quadratic invariant in terms of the quadratic 
fields is universal for all n values. Note, that already for n = 3 one has also higher power group 
invariants. 

The mass-square matrix is evaluated in terms of the derivatives of the invariants with respect 
to the field components and the derivatives of the potential with respect to the invariants. For a 
specific parametrisation one chooses nonzero values for certain Sj variables (si = vi, 7Ti = 0). On this 
background there is no mixing between the s- and the 7r-modes. The genuine expression is identical 
for both, therefore we write the explicit form only for the s-sector: 
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The derivatives of p are the same for all n: 
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From here one can see that the first term on the right hand side of ([6]) contributes to each mode and is 
field independent (equal to 4?7 fc p ). In the indices the quantity appearing after the prime refers to the 
partial derivative of the potential with respect to the corresponding invariant. Then one can write 
the factors k 2 + M 2 appearing in (j^J) as K 2 + p 2 a , K 2 = k 2 + AUk tP , where pi, is the remainder of the 
mass eigenvalue after subtracting from it the field independent piece. In the U(n) x U(n) symmetric 
theory the f-dependent coefficients of the derivatives of Uk start quadratically. The expansion of the 
characteristic polynomial in powers of K 2 then starts as: 

Tl a (K 2 + p 2 a Y a = (K 2 ) N + (K 2 f- l Tip 2 + 1/2 ■ (K 2 ) N - 2 ((Trp 2 ) 2 - Tr(p 2 ) 2 ) + 0(v 6 ), (8) 

which after throwing away terms which contain in the prefactor of the derivatives of Uk at least six 
powers of the f-background, leads on the right hand side of (EJ) to the ratio 



E 



9c 



eigen—a 



k 2 + 4U k , P + p 2 



(9) 



N(K 



2\N-1 



+ (N- l)(K 2 ) N - 2 Trp? + {N — 2)/2 ■ (K 2 ) N -\(Tip 



,2\2 



Tr(^i 



2\2\ 



(K 2 ) N + (K 2 ) N - l Tip 2 + 1/2 • (K 2 ) N ~ 2 ((Trp 2 ) 2 - Ti(p 2 ) 2 ) 
(N = J29a = 2n 2 ). For the derivation of the present form one has to recognize that 

E 9aP 2 a = Tr/i 2 , 9aP 4 a = Tr(/x 2 ) 2 . 
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The naturalness condition for this truncation will be clarified in the section discussing the 3-flavor 
model. 

Assuming further \\fi 2 \\ <C K 2 one arrives at the form which is used to derive the leading order 
perturbative RG-running of the couplings characterising the lowest power operators figuring in the 
power series expansion of Uh [25]: 

k d+l Q k _ / g a g a fj, 2 a g a (n a ) 4 \ . . 

" d (2vr)^_ Q U 2 (^ 4 ) 2 (^ 6 ) 3 ")' { } 

One should note, however, that normally one integrates fl3J) down to k = and in some regions of the 
field variables (especially near the minimum of the potential) the assumed inequality might not hold. 
This circumstance restricts the range where one can rely on solutions of the Wetterich equation in 
terms of a power series. 

In Section 2 the complete expression will be given for n — 2, e.g. without any truncation in 
the characteristic polynomial. In Section 3 we argue that the omitted group invariant makes it 
inconsistent to keep any term of 0(v 6 ) in the expansion of the characteristic plynomial, which yields 
a natural criterium on how and where to truncate. In Section 4 we consider the same truncation in 
case of SU(3) x SU(3) symmetry, where the cubic invariant (det$ + det$^)/2 will be introduced. 
The explicit formulae appearing in these three sections might be known to some workers actively 
pursuing FRG-investigations of the linear sigma model, but to our best knowledge they did not 
appear in previous publications. A short summary and outlook closes this note. 



2 Complete LPA equation for the U{2) x U(2) symmetric 
model 

There is no further invariant beyond those in (J5J) in this case. We choose for the background non-zero 
vq and t>3. The quartic invariant looks then explicitly as 



r = 16 



soSi + iromy + irfsi-(iT l siy ^IQv'vl, k,l = 1,2,3. (12) 



After subsituting the non-zero field values its derivatives have the following expressions: 
dr 2 dr 

— = 32(s s k + n 7i k )s k 32v v 3 , — = 32(s s k + 7i 7i k )ii k 
os dir 

dr 

— = 32[(s s; + ir ni)s + n 2 k si - 2-K k S k -Ki] ->■ 32^3^ 3 , 
dr 

— = 32[(s si + 7r iTi)Tr + s 2 k 7ii - 2n k s k si] -> 0, 
d-Ki 

cP' 

327^0, ^-=32W«, ^-^- = 32^,, 1^3, 



drcodno drcodTTi " dnidrck 

cP 1 cP cP 

32v 2 3 , = e^fe, 1 — l- = 32v 2 6 lj . (13) 



ds ds ' ds dsi ' * dsids k 



The mass matrix is built up readily and one finds directly two degenerate eigenvalues (Z = 1, 2) plus 
the coupled (0, 3) modes both in the scalar and the pseudoscalar sectors. In the matrix eigenvalues 
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one can replace vq and t>3 by the two invariants, but the individual matrix elements will show rather 
non-analytic dependence on r and p (see for instance Eqs.(5.2) and (5.3) of Ref.[13]). Any truncation 
consisting of omitting some of the modes requires particular care in order not to loose the analiticity 
forced upon by the construction of the equation. 

There is no need here to form the full characteristic polynomial of the mass matrix, since after 
some trials one finds that the contribution of the modes grouped into appropriate pairs can be 
expressed in polynomial form through p and r. Actually, for each of these two-dimensional subspaces 
the three terms appearing in dHJ) represent the exact result (there is no 0(v 6 ) contribution). In the 
following formulae we denote the scalar eigenmodes with an extra index s, and the pseudoscalars 
with p. One finds: 

(K 2 + p 2 sll )(K 2 + p 2 pU ) = {K 2 + p 2 s22 ){K 2 + p 2 22 ) = K* + AK 2 + B, 

A = !6pU kjT , B = 64rC/| >T 

[K 2 + p 2 p03+ )(K 2 + ^ 03 _) = K 4 -6ArUl T 

(K 2 + /4 3+ )(^ 2 + /4 3 -) = K* + CK 2 + D, 



C = 8 [2pU k)T + pU kiPP + AfrrU ktTT + 4rU k) 



D = 64 



4t(t- P 2 )(Ulpr-U kjPP U k , TT) 



- U k;T (3rU kiT + 6r 2 U k , TT + AprU k , pr + (3r - 2p 2 )U ktf)p )\ . (14) 

Here the two eigenvalues emerging from the coupled pseudoscalar (0,3) modes are indexed as p p03 ± 
and for the scalar sector the corresponding notation p 2 03 ± is used. 

The pieces in (jT4"l) are substituted into ([2D and one arrives at the full LPA equation. The exact 
invariant representation of the right hand side of this equation in the U(2) X U (2) case is analytic in 
the sense that a power series expansion in integer powers of p and r obviously exists. 

It is clear that keeping only the lightest scalar and three of the pseudscalars, which would reduce 
this model to the usual pion-sigma model, leads to an inadvertent violation of the analiticity of the 
potential in p and r. 

Phenomenological informations suggest that the chiral symmetry breaking pattern does not imply 
any breakdown of the isospin symmetry, e.g. the minimum of the effective potential is expected to be 
located at small values of r (corresponding in our parametrisation to small v 3 ), though its quantum 
fluctuations contribute to the RG-evolution of the potential. This remark suggests the simplifying 
ansatz [131 H3]: 

U(p,r) = U^(p) + U^ 2 \p)r, (15) 

replacing the problem of the determination of a two-variable function by two one-variable ones. 
Technically, one substitutes this ansatz and its appropriate derivatives into fTT4l) and expands the 
expression of the right hand side of ([3D in the variable r up terms linear in r. By comparing the 
coefficients of r° and of r 1 one finds the one variable equations solved in Refs. [T3"l [T4] without 
the need for any extra consideration concerning the cancellation of apparent singular contributions 
unavoidably showing up when this expression is written in terms of vq and v 3 . With a further 
approximation it might be convenient for phenomenological purposes to fix the potential to the 
perturbatively renormalisable form [211 [25] and use Eq. f[TT]) just to trace the scale dependence of 
its couplings[15j. It would be, however interesting to see if there is in the potential any non-trivial 
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structure along the r direction, which can be studied by keeping the full expression on the right hand 
side of ©. 



3 Consistent truncation in the U(3) x U(3) case 

The pattern of chiral symmetry breaking in the three-flavor case is mostly investigated in the v o — v 8 - 
plane. We also used this background for the parametrisation of the quadratic and quartic invariants: 



P = 2(s 2 a + K 



-> 2(v 2 + v 2 ) 



T 



-dij k {siSj + niTTj) + \ -(s s k + n n k ) + {fijk^iSjf 



-vj(v 8 - 2V2v ) 2 . 



(16) 



The derivatives necessary for building up the mass-square matrix are simply obtained for p, and after 
a lengthy, though straightforward calculation one finds the following formulae for the derivatives of 
t: 



— = -v 8 (V2v - f 8 )(2\/2fo - v 8 )Si 8 , 
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Making use of the structural constants of the U(3) group one constructs the mass-square matrix. Its 
scalar and pseudoscalar sectors are independent and of the same structure. The modes k = 1, 7 
are directly found to be eigenmodes, and the modes 0, 8 again form two 2x2 coupled subsystems. 

In the Z7 (3) x U (3) case we could not identify smaller subsets of the factors of the characteristic 
polynomial which could be separately expressed through the invariants, therefore we turn to the 
expansion ([51) . For the terms displayed in this formula we find: 



Tr^ 2 
Tr(/i 2 ) 2 



8pC4,p P + ^pU K r + 32r£4, pr + 0(v 6 ), 

/ 2048 

3SArU k , T U k , pp + Up 2 Ul pp + (512r + —p 2 ) rU 2 KT + 0(v 



(19) 
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In the first line the term 0(v 6 ) refers to the coefficient of Uk, T r, which contains altogether six factors 
of Vo and v$. It receives contribution exclusively from the scalar (0, 8) sector: 

O" + (£f = f T (" + Slg, « " 2 ^»^i/f ) • < 20 > 

This contribution is clearly non-analytic, what is expected to appear just at this level since the 
£7(3) group has a further invariant (~ Tr($T$) 3 ) containing terms with six factors of the fields. In 
the second line of ffT9j) among the omitted terms the lowest number of v j-factors appears in front 
of Uk^prUk^ and Uk,p T Uk,p P - Direct calculation demonstrates that none of these coefficients can be 
expressed through integer powers of the two invariants, i.e. in a form Cip 3 + C2Tp. For this, also the 
group invariant containing six powers of the background should have been included. 

The only consistent way of restricting the dependence of the potential to p and r and obey the 
analiticity is to truncate the expansion of the characteristic polynomial after the three first terms 
and solve the FRG-equation with (jHJ) on its right hand side. 



4 Including the anomaly: the 577(3) x 577(3) case 

An effective linear sigma model which reflects the £7a(1) anomaly is phenomenologically more ap- 
pealing. In phenomenological investigations j2H 125] one does not separate from the matrix field $ the 
zeroth component, which would produce a large number of operators in the FRG program. On the 
basis of processes mediated by instanton configurations one just adds to the U(n) x U (n) invariant 
Lagrangian a term which is less symmetric. The corresponding piece of the Lagrangian is 't Hooft's 
determinant: 5 = (det$ + c.c.)/2. It breaks the U(3) x U(3) symmetry down to 577(3) x SU(3). 
In terms of the parametrisation of the preceding section one has 

5 = A^V^o + v 8 f(V2v - 2v 8 ). (21) 

Below we follow the usual process of taking into account the axial anomaly and just include S into 
the set of variables of the effective potential. We consider U (p, 5, r) and calculate the complementary 
piece to the mass-square matrix of the model with help of the derivatives of S with respect the scalar 
and pseudoscalar fields. These have the following form on the previous background: 
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dsidsj 




d 2 S 
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-2y ~(v - V2v 8 )5 ij , 
2 

-—^(V2v + v 8 )5 i:j , 




dsl 



dn 2 



if i= 1,2,3, 



if i = 4,5,6,7, 



7 



The opposite sign of the second derivatives with respect to the scalar and the pseudoscalar fields 
ensures the cancellation of the contributions to the mass-square matrix which would be linear in the 
background fields. 

The completion of the mass-square matrix involving at least one partial derivative of U k with 
respect to 5 is worked out and the following additions are found to the traces (119]) : 

ATifi 2 = (Jp 2 - r) U m + 2W Mp + 0(v 6 ), 

ATr(/i 2 ) 2 = 32pf/ 2 5 + 1925U kjPp U kj5 - 1024H7 fc)T E7 M + 0(v 6 ). (24) 

One should wonder about the 0(v 5 ) terms of these expressions. In the first line it appears in front 
of Uk,5ri on the second line there is another 0(v 5 ) coefficient in front of U kjPT U k ^. Both are fully 
determined by the (0, 8) scalar sector having the respective coefficient: 

^ / dr 85 dr d5 \ 
\8s 8s ds s ds 8 J ' 

dp dr d 2 5 ( dp dr dp dr \ d 2 5 dp dr d 2 5 

ds ds dsl \ds ds 8 <9s 8 <9s / ds ds s ds 8 ds 8 dsl 

One verifies explicitly using the derivatives listed above that both coefficients vanish! Therefore the 
truncation of the characteristic polynomial will be done again at terms containing no more than four 
factors in the background fields, though the accuracy is still 0(v 6 ). 

It is just the present system where the advantage of the invariant formulation of the RG-equation 
over the use of the component background becomes evident. Although taking appropriate derivatives 
of Uk(vo,vg) one can trace the evolution of the coupling in front of 5, with the vq — v$ component 
parametrisation one can explore only a surface in the three-dimensional p — r — 5 space. On the other 
hand using the invariant formulation, e.g. solving Eqs. (j3J) with fl9]), where the traces are given by 
the sum of (IT31) and (12~4"I) . one can explore the full space. Prospectively, by generalising the proposed 
procedure to finite temperature one can study also the thermal evolution of the expectation value 
of 5, e.g. that of the Ua{X) anomaly. It might give access to the old question whether the chiral 
symmetry or the axial anomaly is restored first [23]. 

The approach described in this note is followed often for 0(n) symmetric theories. There, however, 
a single invariant exists, therefore there is no big difference in using the invariant or the coordinate 
representation of the background. In case of the matrix models with unitary symmetry the formula- 
tion presented here allows a more transparent understanding of the symmetry features of the ground 
state. The choice of appropriate coordinates for the representation of the background is convenient 
to compute the coefficients of the derivatives of the effective potential. However, some cancellations 
(for instance, the vanishing of the coefficients of U k ^ T and U ktPT U ki s) in presence of the determinant 
term appear somewhat miraculous when one uses explicit coordinate representation of the invariants, 
though it might have some simple explanation in terms of group theory. 

Numerical investigations should reveal the possible existence of local minima of the effective po- 
tential which correspond to symmetry breaking patterns beyond the pattern realized in nature. The 
generalisation of the present treatment to the more general case of non-trivial field renormalisation 
is also worth of consideration. 
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